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Abstract: We explore single and multi-loop conformal integrals, such as the ones ap- 
pearing in dual conformal theories in flat space. Using Mellin amplitudes, a large class of 
higher loop integrals can be written as simple integro-differential operators on star inte- 
grals: one-loop n-gon integrals in n dimensions. These are known to be given by volumes 
of hyperbolic simplices. We explicitly compute the five-dimensional pentagon integral in 
full generality using Schlafli's formula. Then, as a first step to understanding higher loops, 
we use spline technology to construct explicitly the 6d hexagon and 8d octagon integrals 
in two-dimensional kinematics. The fully massive hexagon and octagon integrals are then 
related to the double box and triple box integrals respectively. We comment on the classes 
of functions needed to express these integrals in general kinematics, involving elliptic func- 
tions and beyond. 
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1 Introduction 

The scattering amplitudes of M = 4 supersymmetric Yang-Mills theory [1, 2] with SU(N) 
gauge group in the planar (N — > oo) limit are remarkable objects, possessing many non- 
obvious properties. Chief among these are superconformal and dual superconformal sym- 
metries [3, 4] which close onto a larger group of Yangian symmetry [5-8]. Such symmetries 
together with on-shell recursions [9, 10], unitarity-based methods [11-13], the Grassman- 
nian formulation of amplitudes [14-16] and the Wilson loop/scattering amplitude dual- 
ity [17-21] have greatly expanded our understanding of the J\f = 4 theory. These and other 
developments are reviewed in [22]. 

One of the important outcomes of these ideas has been the tremendous progress in our 
knowledge about the structure of multi-loop amplitudes. Although the integrand of the 
theory has been completely constructed [23-26], new mathematical techniques are neces- 
sary to efficiently describe the integrated objects. One advancement along these lines has 
arisen from the study of the six particle Maximally Helicity Violating (MHV) amplitude 
or its remainder function at two loops [27-30]. Dual conformal symmetry, together with 
the proposed duality between scattering amplitudes and Wilson loops, enabled the analytic 



- 1 - 



computation of this particular process, albeit at first only in a very complicated form. How- 
ever, the knowledge about the possible space of functions of the scattering amplitude, which 
in this case were particular types of iterated integrals, and the application of the math- 
ematical tool of symbols particularly suited to such integrals (see in particular [31-33]), 
allowed for the simplification of this remainder function to a very manageable form [34]. 
Symbol technology has seen several other important applications including [35-39]. When 
the number of kinematic variables is sufficiently small (in particular, for several six-particle 
processes), it has even been possible with sufficient effort to obtain analytic formulas for 
certain amplitudes (or Regge limits of amplitudes) [40-44], but our analytic knowledge 
of more complicated integrals appearing in multi-loop SYM scattering amplitudes is still 
quite limited. 

In a seemingly unrelated development it has been realized that the correlation functions 
of Conformal Field Theories (CFT's) in AdS/CFT at strong coupling have properties 
analogous to flat space scattering amplitudes in an auxiliary space called Mellin space, first 
introduced by Mack [45, 46] and further studied in the works [47-53]. An application of this 
formalism in the context of flat space conformal integrals has appeared in [54]. In particular, 
it was shown that a large class of conformal integrals — including those corresponding to 
position space correlation functions in (fr theory, which correspond to various kinds of box 
integrals — have a very simple Mellin representation which can be constructed in terms of 
Feynman rules. Using these, it is straightforward to see that there are simple integro- 
differential relations between various kinds of multi-loop integrals and lower loop ones, 
all the way down to a set of basic building blocks: the one- loop n-gon integrals in n 
dimensions, also known as the n-point star integral. These relations generalize various 
differential relations between integrals of different loop order which have long been very 
useful in the study of scattering amplitudes (see in particular [55-57] for some recent 
examples relevant to SYM theory). 

These results suggest that it is of pressing importance to understand the star integrals 
in detail (a close relative of our star integrals, with massless external legs but massive 
propagators, has been studied and evaluated explicitly in several cases in [58]). In this 
note we take some modest steps in this direction. Firstly, it has been realized that such 
integrals compute volumes of simplices in hyperbolic space [59-62] (a different relation 
between amplitudes and volumes has been explored in [16, 63]). We can therefore use 
Schlafli's formula, which determines the differential of the volume of an (n — l)-simplex in 
terms of the volumes of (n — 3)-simplices (a motivic version of Schlafli's formula [64] has 
been similarly applied to compute symbols of star integrals in [65]). As one application, 
we integrate the formula explicitly to find the d = 5 pentagon integral. The result is 
remarkably simple, being simply a sum of logarithms with unit coefficients. The d = 6 
hexagon and d = 8 octagon are addressed next. In these cases finding the full answer 
appears much more difficult (some special cases of the d = 6 hexagon have been explicitly 
evaluated in [56, 66-68]) and we will content ourselves with finding analytic results when 
the external kinematics are restricted to two dimensions. We apply the results of the 
recently developed spline technology for loop integrals [62], which tells us that in such 
kinematics, these integral can be written out as sums of box integrals with determined 
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coefficients. 

The fully massive d = 6 hexagon (d = 8 octagon) integral plays a role in determining 
the fully massive double (triple) box integrals in four dimensions. The relation of the d = 6 
hexagon to the double box has been worked out in [54]. In this note we do the same for 
the triple box and the octagon. We find the former is given as a double integral of the 
latter. Crucially, the hexagon and octagon integrals being integrated over are ratios of 
polylogarithm functions divided by certain square roots. We argue this has implications 
for the class of functions in terms of which higher loop integrals can be expressed. 

The paper is organized in the following way. In Section 2 we review general ideas 
about Mellin amplitudes and the consequences of the existence of Feynman rules for Mellin 
amplitudes. Namely, we discuss the connections of multi-loop Feynman amplitudes with 
products of Mellin amplitudes, and its implications for the position space results. We 
stress that from this it is clear that we need to have a better understanding of ra-gons in n 
dimensions or the "star" integrals to understand the fully massive loop integrals of J\f = 4 
SYM. In Section 3 we discuss these star diagrams in more detail reviewing some known 
results as well as presenting some new analytic results for pentagons in five dimensions. 
For more complicated diagrams like the d = 6 hexagon and the d = 8 octagon it is very 
difficult to get explicit results for general kinematics. So, in Section 4 we extensively discuss 
the analytic results for 2n-gons using a restrictive kinematic localized in two dimensions. 
We do this by using the technology of splines to simplify such computations and present 
explicit results for two examples, the d = 6 hexagon and the d = 8 octagon. In Section 5 
we determine the representation of the triple box integral as a double integral of the d = 8 
octagon. Both in this case and for the double box, the integrand has a square root in 
the denominator which we know explicitly. We study various kinematic limits which tell 
us whether or not one should expect to see elliptic, or even more complicated, functions 
rather than the generalized polylogarithms which are much more familiar in multi-loop 
computations. Our results agree with the analysis of [69]. Some details about our results 
for the d = 6 hexagon and the d = 8 octagons in 2d kinematics are collected in the 
Appendices, and a Mathematica notebook with the full expressions is available in the 
online version of this note. 

2 Mellin amplitudes refresher 
2.1 The Mellin amplitude 

The multi-dimensional Mellin transform formalism was introduced in the work of Mack [45, 
46] and quickly applied to both AdS /CFT [48-53] and flat space calculations [54, 70]. The 
Mellin transform can be applied to any conformally invariant function of several points Xj, 
with given conformal weights Aj. This could be a conformally invariant correlation function 
or a conformally invariant integral (in applications to SYM theory scattering amplitudes, 
these are usually called dual conformal as a reminder that the relevant conformal symmetry 
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is that in momentum space, rather than position space). For instance, we can write 




where Xy = (xj — Xj) 2 and the 6ij parameters satisfy the constraints 

^<% = 0, 5 U = -Ai. (2.2) 

The function M(5ij) is usually called the Mellin transform of {4>Ai • ■ ■ <M„ (x n ))- After 
solving the constraints, the integral becomes an ordinary multi- variable Mellin transform 
in terms of n(n — 3)/2 independent variables. The integration is over a set of complex 
variables Cj, each running from — ioo to +ioo along an appropriate contour. The con- 
straints (2.2) guarantee that the variables Xij in the integrand combine into cross-ratios, 
thereby imposing conformality. It is important to note that the constraints can formally 
be solved by introducing a set of Mellin momenta ki, satisfying momentum conservation, 
J2i ki = 0, such that 

$ij = ki ■ kj, k\ = A.;. (2-3) 

This parameterization provides some intuition for the 5ij parameters. In fact, in practice 
it is convenient to work with Mandelstam type variables, s^...* = —(k^ + . . . + ki p ) 2 , e.g. 
si2 = ~{ki + k 2 ) 2 = Ai + A 2 - 2«5i2. 

2.2 Feynman rules and convolutions 

In [54] a subset of us found that Mellin transforms of the kind of (dual conformally in- 
variant) integrals that appear in SYM theory scattering amplitude computations have an 
extremely simple form. Consider for example a momentum space diagram whose position 
space dual is the same as a position space correlation function in </> 4 theory (three exam- 
ples are shown in Figure 1, with the dual graphs shown in blue). The Mellin amplitude is 
obtained from the dual graph by the simple rules: 

• To each external leg associate a Mellin momentum ki such that kf = —1. 

• Momentum flows through the diagram being conserved at each vertex. 

• To each internal leg with momentum k associate a propagator l/(k 2 + 1). 

In other words, the Mellin amplitude looks just like a momentum space amplitude for 
massive <p A theory, with m 2 = 1. This 1 is nothing but the canonical dimension of <fi, 
A = (d-2)/2 = 1. 
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Figure 1. The one-, two- and three-loop ladder diagrams (black) and their corresponding dual 
tree diagrams (blue). The external faces of the former, or equivalently the external vertices of the 
latter, are labeled xi,X2, . . . clockwise starting from x\ as indicated. 



According to these rules we have, for example, the following very simple results for the 
Mellin amplitudes of the box, double box, and triple box integrals shown in Figure 1: 



M = 1, 
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The Feynman-like rules nicely express Mellin amplitudes as products of simple factors. 
We can use this to our advantage since a product in Mellin space maps back into position 
space as a convolution of the individual position space expressions. That is, suppose we 
have two functions f(x),g(x) with Mellin transforms M* (s), M 9 (s), 



M f (t 



+oo 



dx 



x 



x s f(x), M 9 {s) 



dx 



x s g(x). 



(2.7) 



Then the position space representation for the product M$ (s)M 9 (s) is 

h(x) = <b — Mf(s)M 9 {s)x- s 1 ds f+C ° dy 
J 2m 

— f(y)g(x/y)- 
o v 



ds 
2tH 



y s f(y)M 9 (s)x- 
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Accordingly, we can split the computation of higher-loop integrals into two steps: first we 
compute the position space expression corresponding to the Mellin transform, which is just 
a product of propagators; and the second, more difficult step is to evaluate the position 
space expression of the product of T functions appearing in (2.1). But the latter is nothing 
but the same as computing a diagram whose Mellin amplitude is M = 1, which corresponds 
to the n-legged star graph, examples of which are shown in Figure 2. 

In SYM theory amplitude calculations we are also often interested in diagrams with 
various numerator factors. These can be translated into Mellin space as differential oper- 
ators acting on the Mellin amplitude. Therefore we expect that a large class of integrals 
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Figure 2. The 'star' graphs for n = 4,6, 8, in blue, correspond to the one-loop box, hexagon, and 
octagon integrals in d — 4, 6, 8 respectively. These are the basic building blocks for many integrals 
relevant to multi-loop scattering amplitudes in SYM theory since each one is simply M = 1 in 
Mellin space. 



which appear in SYM theory scattering amplitude computations, to all loop order, can be 
expressed as integro-differential operators acting on just one class of elementary object: 
the n-point star integral in position space <j) n theory, or equivalently the one-loop n-gon 
Feynman integral in n dimensions. This makes it clear that studying these objects is an 
important first step in understanding the analytic structure of a large class of multi-loop 
integrals. 



3 Star integrals 

It is convenient to use the embedding formalism [71, 72]. This amounts in practice to 
defining d + 2-dimensional null vectors P M to describe d-dimensional coordinate vectors 
x 11 , via 

P M = (P+,P-,P») = (l,x 2 ,x^). (3.1) 

It is easy to check then that P- L j = —2Pi ■ Pj = (xi — Xj) 2 = xfj. 
The n-gon star integrals are defined by 

/w = r ddx n 1 - f ddQ fr 1 (32) 

J i7r d/2 11 {xi _ x) 2 J in d/2 11 { _ 2Pl . Q) ■ ^ ■ > 

They are simply related to volumes y( n_1 ) of ideal hyperbolic (n — l)-simplices [59-62] 
according to 



y(n-i) _ y_] detPjjl ^ (n) 



»r«) (3 ' 3) 



Let us now consider the first few cases. 



Triangle 

It is straightforward to do the integral directly in this case, and one finds 

7 (3) = F (2) 3 (34) 
V -f 12 Pl3 P23 

Using formula (3.3) above this gives = tt, which is indeed correct: the area of a 
hyperbolic ideal triangle is precisely equal to tt. 
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Box 

The simplest non-trivial star integral is the first one in Figure 2, corresponding to the 
four- dimensional box function. The result for this well-known integral is given by 



Li 2 (x + /x_) - Li 2 + Li 2 -(x + ^x 



in terms of 



-7"^ Jdetx^ 



1 



(3.5) 



x± = - ( 1 + ui - u 2 ± \l 1 - 2u\ + u\ — 2u 2 - 1u\U2 + u\ J (3.6) 



and the two cross-ratios 
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The numerator in (3.5) is nothing but the Bloch-Wigner function (see e.g. [61]), which 
indeed is known to compute the volume of an ideal hyperbolic tetrahedron. 

Pentagon 

The next-simplest case, not shown in Figure 2, is the one-loop pentagon integral in five 
dimensions, which as far as we are aware has not been explicitly evaluated in the literature 
(the one- loop pentagon integral in four dimensions has been evaluated in [73]). Surprisingly, 
we find that it takes a very simple form. 

The pentagon integral corresponds to the volume of a hyperbolic 4-simplex. Such a 
volume depends on five cross-ratios, which in turn are built out of the five coordinates a^. 
Let us take concretely 

f*14 -P23 -P25 -P34 -Pl3 -Pl5 -P24 Pl2 P35 ro x 

U1 = ^W U2 = ^W U3 = iWV U4 = ^' U5 = PnPV 5 - (3 - 8) 
To obtain an expression which only depends on cross-ratios we consider the rescaled integral 



/ (5) = \^P^APV 5 I {5) . (3.9) 

The computation of the volume is most straightforwardly done using Schlafli's formula. 
The formula relates the differential of a hyperbolic simplex in terms of its co-dimension 2 
simplicial faces and associated angle differentials — since each co-dimension 2 face is defined 
by the intersection of two hyperplanes (which lie along co-dimension 1 faces), there is 
therefore an associated angle. This angle can be represented in terms of the vectors normal 
to said hyperplanes. 

More concretely, if we have a simplex whose vertex representation is given by the Pi 
vectors, its hyperplane representation is given in terms of vectors Wi which are normal to 
these hyperplanes. In particular, Wi ■ Pj = In terms of these we can write Schlafli's 
formula as 



-1 » m ... (Wi-Wj + Jw-Wjy-wfw. 



2*(*-l)^ ( fe 2 ) I Wi ■ Wj - J {Wi ■ Wj) 2 - Wjw] 
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(a) 

where V\_ 2 ^ corresponds to the volume of the d — 2 simplex spanned by all the vectors 
except for the pair Pi, Pj. 

This formula is particularly simple in the case k = 4. In this case the V(fc-2) become 
volumes of ideal hyperbolic triangles. But this is simply ir\ The integration of Schlafh's 
formula depends on the kinematic region under consideration. We work in the Euclidean 
region where all (xi — Xj) 2 are positive, and if we define 



A® 



1 



det Pi 



2 P13 Pl4 P24P25P35 

1 - [wi(l - ii 3 (l + U4) + u 2 u\) + cyclic] - uiu 2 u 3 U4U 5 , 



(3.11) 



then for A^ 5 ) < we have 



V, 



(4) 



£ ( 

l<i<j<n 



■lY+i lOE 



Wi ■ Wj - J(Wi ■ w 3 ) 2 - w 2 w 2 



Wi ■ Wj + J(Wi ■ w 3 ) 2 - w 2 w 2 



(3.12) 



and using (3.3) this gives 



f(5) 
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7T2 



2v/-A( 5 ) 




-l) i+J 'log 



Wi ■ Wj - ^J{Wi ■ w 3 ) 2 - wfw- 



Wi ■ Wj + J(Wi ■ Wj) 2 - W 2 Wf 



The pentagon has a cyclic permutation symmetry under the action of g : Ui 
can then finally write the remarkably simple and manifestly symmetric form: 




/(5) 
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7T2 



--{l + g + g 2 + g 3 +/)W 



V-A( 5 ) 
V-A( 5 ) 




s + 



V-A( 5 )' 
V-A( 5 ) 



with 



2V-A( 5 ) 

(1 - U 2 )(l - U 5 ) - Ul(2 — U S — W4 — - U 2 ""4 + «lM3«4) 

r = 2 ' 

(1 - u 5 )(l - u 2 u 5 ) - ui (1 + «5 — 2 n 3 u 5 + u 4 + u 2 ^4^5 + ^1^4) 



>(3.14) 

(3.15) 
(3.16) 



2a/Ui~U5 

Hexagon and beyond 

Using Schlafh's formula (see [64, 65] for further details), one can easily express the differ- 
ential (or, if one likes, the symbol) of the ra-dimensional re-gon integral as a sum of certain 
n — 2-dimensional n — 2-gons. However, it is in general a difficult task to integrate this for- 
mula analytically. The structure of the differential equation makes it clear however that it 
can always be expressed in terms of generalized polylogarithm functions [64]; in particular, 

• i"( 2n ) can be expressed in terms of functions of transcendentality degree n, 

• and jC 2 ^ 1 ) can be expressed in terms of functions of degree n — 1. 
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One way to understand the apparent inconsistency of the transcendentality counting in the 
two cases is that the odd-dimensional integrals always contain an overall factor of 7r 3 / 2 , 
as we saw explicitly for the pentagon in (3.14). Taking this factor into account, the Tri- 
dimensional m-gon integral always has total degree m/2. We remind the reader that all 
generalized polylogarithms of degree less than 4 can be expressed in terms of the classical 
polylogarithms Li m , so non-classical polylogarithms first appear in the d = 8 octagon 
integral (for general kinematics). 

We turn now to the d = 6 hexagon integral, which has received attention in the 
literature [56, 65-68] in part due to its interesting relationships (via differential equations) 
to other integrals relevant to SYM theory scattering amplitudes [56] . However it remains an 
interesting outstanding problem to fully evaluate the d = 6 hexagon in general kinematics, 
where the integral depends on 9 independent cross-ratios (we present a choice of cross- 
ratios in Appendix A). To date the closest we have to this is the analytic formula for the 
special case of the "three-mass easy" hexagon [68] (an expression for its symbol was given 
in [65]). In this case three of the nine cross-ratios are set to zero. The formula presented 
in [68] therefore computes the d = 6 hexagon on a six-dimensional subspace of the full 
nine-dimensional cross-ratio space. 

Motivated by the desire to simplify the evaluation of otherwise difficult integrals, and 
by the vast body of recent work on SYM theory amplitudes in two-dimensional kinematics 
(see for example [57, 74-77]), in this paper we therefore carry out explicit computations of 
the d = 6 hexagon and the d = 8 octagon in 2d kinematics. Here, due to Gram determinant 
constraints, the nine cross-ratios for the hexagon (and the twenty cross-ratios for the general 
octagon) are constrained to take values in a six-dimensional (ten-dimensional) subspace of 
the full parameter space. We present explicit parameterizations of the cross-ratios in terms 
of six (ten) free variables in Appendices A and B. Our result for the d = 6 hexagon in 
2d kinematics is in a sense complementary to that of [68] since the two six-dimensional 
subspaces are disjoint inside the full nine-dimensional parameter space of the generic d = 6 
hexagon. 



4 2n-gon loop integrals in 2d kinematics 
4.1 Setup: splines 

In this section we evaluate 2ri-gon loop integrals in two-dimensional kinematics. To do this 
we shall use the methods developed recently in [62] based on spline technology, which the 
reader should consult for further details. With these, it can be shown that the one-loop 
star integral (3.2) can be written in the form 



j(n) 

where the spline is defined by 



2 / e x2 T(X;{P t }) (4.1) 
Jm D 



T(X;{P l })= / HdnSW(X-J2tiPi)- (4-2) 

"'° i=l i=l 
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This expression follows by noticing that the spline is the Laplace transform of the integrand. 
Here we are interested in 2d kinematics, so we set D = d+2 = 4. We shall also only consider 
even-dimension integrals and therefore set n — > 2n. The computation of the spline depends 
on the various linear relationships between the -Pi's. Here we shall assume that the vectors 
are generic, i.e. that every set of four vectors spans M 4 . 

Under these conditions the spline can be written as a sum of terms, each corresponding 
to a particular linearly independent set of vectors. Not all such sets need be considered 
though. It is sufficient to take the set B of so-called unbroken basis, which for generic 
kinematics amounts to the set of basis which include the vector P\. To each such basis, b, 
there corresponds a piece in the spline, which is therefore made up of N = (n — 1)1/ (n— 4)!3! 
terms. Each term is labeled by its unbroken basis, b, and the coefficients can also be easily 
computed. In this manner we find 

. xm(X) 

'^•"-& ie5<w»>.P»> >/*«<n- (4 - 3) 

Some explanations are in order. Firstly, denotes the ith vector not in the basis b. 
Secondly the vectors wf^ are defined by 

wj b) .P j = 6 ij , VPjGb. (4.4) 

We can think of b itself as a matrix whose columns are the vectors Pi G b. This allows 
us to compute the determinant. Finally, X(b){X) is the characteristic function of the cone 
spanned by the vectors in b, which can be written as 



X(b) (X) = l[e(W l {b) -X). (4.5) 

1=1 



To proceed we must evaluate the Gaussian- type integral in (4.1). We could evaluate it 
directly, since the spline is homogeneous in \X\ = y/^C 1 . This would give us a sum of 
integrals of X polynomials over AdS tetrahedra. However, instead of doing this we can use 
the presence of the exponential to integrate by parts the terms of the form W ■ X. At the 
end of this procedure, there are no such factors left, but there are however several types of 
terms, depending on how many times we differentiate the characteristic functions X(b)(X)- 
In particular, one set of terms does not involve derivatives of at all: 



(4.6) 



( (b)\ n ~ 

lin) = (n-4)!! v ) r &x , X {b) {X) 

2t ~ 2 h^B lir=i 4 w i ] ■ ^ M4 Vdet b T b 



This is interesting, as the integrals above are nothing but box integrals, with four external 
legs Pi corresponding to the elements in the basis b. Accordingly, the kind of terms above 
are simply a sum of box integrals, namely dilogarithms. In contrast, the . . . represent terms 
which have an even number of derivatives of X(b)(X)- We have explicitly checked that all 
such terms cancel between themselves for n = 3,4. To understand why, notice that those 
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terms involve for example integrals over lines in AdS, which leads to single logarithmic 
terms. In order to have an expression of uniform transcendentality, it must be that these 
terms actually add up to zero. 

4.2 Applications: hexagon, octagon, and beyond 

To see in detail how we can perform the computation of these coefficients, let us set n = 3 
and consider the particular basis made up of elements Pi, P2, P3, P4. We then have 

C M pN pP pQ 
^(1234), M _ t NPQ r 2 r 3 r 4 



eABCDPfPiP^P? 



W^ 1234 -^ zMNPp p p pApBpC 

1 / _ d ABC P2,MP3,NP4,pP2 P3 4 



234) • P 5 ) (Wi 12 ^ ■ P 6 ) C P 5 ,mP2,nP 3 ,pP4, Q P 6 A P2 B P 3 C P? 



(4.7) 



with " b N the totally antisymmetric product of N delta functions. It is important to 
notice that this expression, when multiplied by the inverse of V det b T b, will have total 
homogeneity —1 in each of the vectors P{. Although we have focused on a particular term, 
this is a generic feature. It guarantees that, if we multiply l2 n by -P14-P25-P36] each term 
in the sum is separately conformally invariant, and can hence be written in terms of the 
nine cross-ratios of a conformal six point function (though we must keep in mind the result 
is only valid for 2d kinematics, which imposes non- linear relations on these cross-ratios). 
We give a choice for these in Appendix A, together with a 2d kinematics parameterization 
for them in terms of 6 independent variables \f ) i = 1> 2, 3. In terms of the latter, we can 
write the contribution of the particular basis (1234) to Iq as 



(2n-4)H XiX 2 xf 



2"~ 2 ((xi - X2 ) xt + {xi + 1) X2xt) {-xf + xt + X! (xt + 1)) 

X B + .. 



(xr + i) (xt-xt) (xi + 1) 2 



D 




(4.8) 

Overall, there are a total often such terms. The total result is too cumbersome to reproduce 
here, but in the online version of this note we include a Mathematica notebook with the 
full result. 

The computation of the d = 8 octagon integral in 2d kinematics is entirely analogous to 
what we have just done. There are now a total of 35 terms in the spline, each corresponding 
to a box integral with a certain coefficient. The d = 8 octagon depends on 20 cross-ratios 
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which in 2d kinematics can be parameterized in terms of 10 independent parameters. The 
details of this kinematics have been included in Appendix B. The full expression for /( 8 ) 
for the d = 8 octagon have been included in the attached Mathematica file since it is very 
lengthy. 

It is straightforward to consider generalizations of the results above and consider 2n- 
dimensional integrals in 2m kinematics, for n > m + 1. Under such circumstances one 
finds the 2n-dimensional integral decomposes (for generic 2m-dimensional kinematics) into 
a sum of (2n — l)!/(2n — 2m)!(2m — 1)! 2m-integrals with well defined coefficients. For 
instance, the general even-dimensional integral in 4d kinematics is given by 

,(6)^ 2n " 6 



(2n) _ (2n-6)!! ^ gTj f x , X{b) (X) 

For the d = 8 octagon the number of unbroken basis made up of six vectors is 21 and 
accordingly the d = 8 octagon is a sum of 21 d = 6 hexagon integrals. 

5 Elliptic functions and beyond 
5.1 The double box 

One of the motivations for this work was to make an attempt to begin exploring integrals 
which evaluate to functions outside the class of generalized poly logarithm functions. Ellip- 
tic functions of this type have been encountered before in explicit QCD computations [78], 
and have been argued to appear in SYM theory as well starting with a double box integral 
contribution to the 2-loop 10-point N 3 MHV amplitude [69]. 

Using the convolution tricks explained in Section 2.2, it was shown in [54] that the 3 
to 3 exchange diagram in position space of (j) 4 theory, which is the same as the double box 
Feynman integral, can be expressed as a one-fold integral of the 6-point star (the d = 6 
hexagon integral): 

I 3>3 (ui,...,u 8 ,u 9 ) = ^/ (6) (ui,...,U8,U 9 ). (5.1) 

Jug u 8 

with /( 6 ) = x^x^x^gl^ and the double box integral, 

j _ f d 4 x a d 4 Xfc xf 4 X2 5 ^3 6 (&2) 

' / ( <l/TJ-2\ ,y> 2 ,y. 2 ,-y> 2 ,y.2 ™2 ,-y> 2 ,>i 2 

J \ U[ ) x la x 2a x 3a x 4fe x 5fe x 66 x a6 

Thanks to our results in Section 4 we are now in possession of a simple formula giving 
the d = 6 hexagon in 2d kinematics. One therefore may hope that this should suffice for 
determining the double box in the same kinematical regime. However the formula above 
demands that the integration is done keeping all cross-ratios fixed except one, and it is 
easy to check that this is impossible in 2d kinematics, since the number of independent 
cross-ratios in this case is reduced because of Gram determinant identities. It is somewhat 
unfortunate that in order to recover a lower-dimensional kinematics result we have to take 
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a detour through the full, generic result. Similar remarks hold for higher loop integrals: 
although we are only interested in Ad kinematics at the end of the day, our convolution 
formulae nevertheless require a detour through a higher-dimensional regime. 

The symbol of the fully general d = 6 hexagon is known [64, 65], but it is rather 
complicated, and integrating it in general remains an interesting open problem. Since 
obtaining the full result seems to be currently out of reach, what can we say about it? 
Well, firstly we know what form the final expression has to take. We know that the d = 6 
hexagon integral is related to the volume of a 5-simplex living in an AdS§ submanifold of 
AdS-j. Denoting this volume by V5, we have from formula (3.3) (and neglecting numerical 
factors) : 

I^)( Xl )^—^=. (5.3) 
' det x\- 

Schlafli's formula tells us that the differential volume of the 5-simplex is fixed entirely in 
terms of that of the 3-simplex, and from this we know that the result will take the form 

/^)^^%i±^, (5.4) 
det x?- 

where in the numerator of course Li3 is shorthand for various terms of the correct tran- 
scendent ality, such as L12Q logQ, log() logQ log() and C(3), with complicated functions of 
the cross-ratios as arguments. 

Our expression for the double box integral then becomes 

h,3 = / — t=^f • (5.5) 

Ju 8 u s VA(«) 

with A^ 6 ) = , , 2 e 2 X % . 2 . In general, A( 6 ) is a third-order polynomial in u%, 

^(6) _ [4 Uin2n5li6U7Ugn 3 _|_ l ower -order terms in u%\ . (5-6) 

Therefore, if any three cross-ratios are set to zero (and both U3 and U4 must be included 
in the three), then the determinant necessarily reduces to a second-order polynomial in u%. 
This is important, since the order of the polynomial determines whether we should expect 
elliptic functions to appear in the final expression for the double box after integrating (5.5). 
Indeed, if we get rid of the polylogarithms for a second, the integral 



/ 



dU8 (5-7) 



u s ^(u s - a)(u 8 - b)(u 8 - c) 



leads to elliptic functions for generic a, b, a If any pair of roots degenerates, or if the 
polynomial becomes second order instead of cubic, we would obtain logarithms instead. 
Because of this, it seems almost certain that the final integrated expression for the double 
box will contain elliptic functions, in general kinematics. 
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Let us look at a particular limit of the general kinematics where we actually expect to 
start seeing the elliptic functions in the final result. For the d = 6 hexagon the "minimal 
massive" case where we go beyond polylogarithms would be the case of 4 massive legs. Say 
we have = and x| 4 = 0. In this case we have U3 = = 0, and as argued above this is 
the largest number of vanishing cross-ratios we can have while staying within the realm of 
elliptic functions. This configuration is exactly the case appropriate to the 10-point double 
box integral shown in Figure 6 of [69]. Now if we further set the other cross-ratios (apart 
from ug) to some constant, generic values, then (5.7) certainly gives an elliptic function, so 
we would expect the same to be true for the double box in (5.5). However any other case 
with a smaller number of massive legs only gives polylogarithms, never elliptic functions, 
because for such cases the polynomial inside the square root degenerates from cubic to at 
most quadratic order. It might be interesting (and certainly easier) to derive the hexagon 
integral with u$ arbitrary, M3 = 114 = and all other cross-ratios set to some carefully 
chosen kinematic values. This would be sufficient to plug into formula (5.5) and check if 
elliptic functions actually do occur there. 

5.2 The triple box 

Let us now consider the triple box integral. In the dual position space this looks like a 
tree-level diagram involving 8 particles and two internal propagators (shown in Figure 1). 
Accordingly we expect it to be given by a two-fold integral of the 8-point star integral. 
This is what we shall proceed to show just now. 

To begin with, we need a basis of cross-ratios which can describe a conformally invariant 
function of 8 points. For fully generic kinematics (in general dimension) we expect 8x5/2 = 
20 independent cross-ratios. We list a choice of such cross-ratios in Appendix B.l. Next, 
we consider the Mellin representation of the triple box. According to the rules we set out 
in Section 2 it is the product of two propagators. Once the constraints (2.2) are solved, we 
get an ordinary multi-dimensional transform in terms of the 20 independent cross-ratios. 
In this way we find 

r +ioo ( 20 1 \ , 8 

- L (n £ <j (2+2ci2 )( 2+2C9 ) * nw m 

with 

j _ 1 f dxgdxbdxg x1 5 X2qx1 7 x1 8 

' "~ 1 / ( ITT^S ry-» 2 ,-y> 2 ,-y> 2 ,-y> 2 ™ 2 /y»2 ™2 ,-y»2 ,yi 2 ,-y> 2 

\ ul ) ■ L la x 2a x 3a x 4b X 5b X 6c x 7c x 8c- l ab X bc 

For the reader's benefit we provide in Appendix B.2 an explicit formula for the product 
of 28 T-functions written out in terms of 20 independent variables q. 

We don't need to display all those details here since we already know that the position 
space expression corresponding to the product of gamma functions is nothing but the d = 8 
octagon integral. We have therefore only to compute the (much simpler) position space 
expressions corresponding to the propagator factors. For instance, 

•+*~ dcis = 8(1 -ma) 

2vri 2(l+ci 2 ) 2ui2 { ' 1 
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with Q(x) = 1 for x > 0, and zero otherwise. In this manner we conclude that 

^3.2,3(^1, •■• ,u 2 o) = / du' g / du' 12 / (8) (ni, . . . ,u' 9 , . . . ,u' 12 , . . . ,u 20 ).(5.11) 

U9 UU J U g J Ul2 

with /( 8 ) = x 2 l5 x 2 26 x 2 37 xl 8 I^ . Of course, this equation can be turned around to write a 
(very simple) differential equation expressing the octagon as a second derivative of the 
triple box. 

Finally let us remark on the d = 8 octagon integral. It is again given by the volume 
of a hyperbolic simplex, and accordingly we have something of the schmatic form 

with A^ 8 ) now given by y-^ — 2 — ^ 2 \2 ■ Of course we emphasize that the numerator will be 
a linear combination of (generalized) polylogarithm functions of degree four, including not 
just L14Q but also for example Lia^Qj Li2()Li2(), etc. If we evaluate the determinant, we 
find 

U13 Uu Uis U19 

where the . . . stands for terms of lower degree in ug or u\2- For general kinematics the first 
integration with respect to ug will make elliptic functions appear, while it is reasonable to 
expect that (again, in general kinematics) the second integration will lead to an even more 
complicated class of functions, beginning as we see at three loops. 
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A Details on the double box computation 

A list of 9 multiplicatively independent cross-ratios required to describe conformally in- 
variant functions of six point is given by the following set: 



,-y> 2 ^-.2 ™ 2 ,-y> 2 o^2 ™ 2 /-y»2 ™ 2 

_ x 14 x 23 _ x 15 x 24 _ x 16 x 25 _ x 25 x 34 

U\ 9 9? ^2 9 95 "3 9 9' ^4 9 9' 

ry" ry " ry " ry" ry " 171^ ry" ry" 

x 13 x 24 x 14 x 25 x 15 x 26 x 24 x 35 



ry" ry" ry" ry" ry" ry" ry" ry" ry" ry*-* 

„ _ x 26 x 35 „. _ x 12 x 36 „, _ x 36 x 45 _ x 13 x 46 _ x 14 x 56 /a i\ 

U5 — —r g - , U 6 — — 2 g - , U7 — — 2 n~, Ug — —~ 2~, Ug — — = „— . (A.1J 

ry" ry" ry " ry" iy " ry " ry " ry " ry" ry" 

x 25 x 36 x 13 x 26 x 35 x 46 x 14 x 36 x 15 x 46 



In order to carry out the computation for the d = 6 hexagon in 2d kinematics we can 
first restrict the general kinematics of (A.l) to a four-dimensional sub-space parameterized 
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by 12 momentum twistors [79]. Subsequently, we can further restrict the 4d momentum 
twistors to a subspace of 2d kinematics which can be very simply parameterized using 6 
independent cross-ratios, as a generalization of the parameterization used in [57, 75, 77]: 
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(A.2) 



/ 



In terms of these variables one may compute the xfj = det(^2i-i-^2i-^2j-l-^2j)) so that the 
9 cross-ratios (A.l) are given by 

Xl" (X3~ + 1) (xi + - X3 + ) 



Ui 



U2 



U i 



"5 



u 7 



us 



Itg 



(xi +i)xa {xi + -x2 + y 

(xT + 1) (X2 + + 1) 

(X3- + 1)(X3 + + 1)' 
(X2~ ~ X3~)X2 + (X3 + + 1) 
(X2~ - (X2 + + 1)X3 + ' 

X3+ + 1 

Xi + Xi" +Xi" +Xi + + 1' 

(X2~ - Xl~) (Xl + + 1)X3 + 

X2"Xl + (X3 + + 1) 
X2~ (Xl~ - X3~) Xl + (X2 + - X3 + ) 



(X2 -Xl )X3 (Xl + 
X2 _ Xl + 



X2"^JX3" 



(X2- + 1)(X1 + + 1)' 
(X2~ + 1)X3~ (Xl + - 



X2 



X2" (X3~ + l)Xl + 
X3~ + l 
X2 + X2" +X2" +X2+ + 1' 



(A.3) 
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B Details on the triple box computation 



B.l Cross-ratios for eight-point functions 



A list of 20 multiplicatively independent cross-ratios required to describe conformally in- 
variant functions of eight points is given by the following set: 



Ul = 


X 15 X 24 

2 2 ' 
x 14 x 25 


U 2 


x 16 x 25 

2 2 ' 
x 15 x 26 


^3 = 


x 17 x 26 

2 2 ' 
x 16 x 27 


ti4 = 


x 26 x 35 

2 2 ' 
x 25 x 36 


u 5 = 


x 27 x 36 

2 2 ' 
x 26 x 37 


Uq -- 


,-y»2 ™2 

_ x 28 x 37 
2 2 ' 
x 27 x 38 


U 7 = 


/T^2 ,-y,2 

x 37 x 46 

2 2 ' 
x 36 x 47 


U 8 = 


™2 ™ 2 

x 38 x 47 

2 2 ' 
x 37 x 48 


Ug = 


x 13 x 48 

2 2 ' 
X 14 X 38 


UlO 


/T^2 ,-y,2 

_ x 48 x 57 
2 2 ' 
x 47 x 58 


U U 


™,2 ™2 
_ x 14 x 58 
2 2 ' 
X 15 X 48 


Ul2 


o^2 ™2 
_ x 15 x 68 

x 16 x 58 


U13 = 


x 13 x 28 

2 2 ' 
x 12 x 38 


U U 


,71 2 ,-y> 2 

_ x 13 x 24 
2 2 ' 
x 14 x 23 


Ul5 


_ x 24 x 35 
2 2 ' 
X 25 X 34 


U±6 


™2 ™2 
_ x 35 x 46 

x 36 x 45 


Un = 


ry 2 ryi2 

x 46 x 57 


U18 


/T-<2 '2, 
_ x 57 x 68 




,-y>2 ™2 

_ x 17 x 68 


U20 


o^2 ™2 
_ x 17 x 28 


2 2 ' 
x 47 x 56 


2 2 ' 
x 58 x 67 


2 2 ' 
x 16 x 78 


x 18 x 27 



(B.l) 



As in the previous section for the double box computation in 2d kinematics we can use 
the momentum twistor parameterization of the above cross-ratios in terms of 16 momen- 
tum twistors in a four-dimensional subspace, which are again expressed in a 2d subspace 
parameterized by 10 cross-ratios. The momentum twistor representation is given by, 
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(3.2) 
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In terms of (B.2) the 20 cross-ratios then take the values 
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(xr + i)(x2 + + i) (xi~ -X2~)x2 + (x3 + + i) 

{Xl~~X2~)X2 + ' (Xl~ - X3~) (X2 + + 1)X3 + ' 

(Xl~ - X3~) X3 + (X4 + + 1) X2~ (Xl~ - X3 _ ) (Xl + - X2 + ) X3~ 

u 3 = 7 Z ZZ7 Z ; T\ T' u 4 ~ 



(Xl -X4 )(X3 + + 1)X4 + ' ' (Xl -X2 )X3 X2 + (Xl + - X3 + ) ' 
X3~ (Xl _ - X4~) (Xl + - X3 + ) X4 + X4~ (Xl~ - X5~) (Xl + - X4 + ) X5~ 



^5 = 7 Z Z\ Z T7 Z ZZ> n 6 



(Xl -X3 )X4 X3 + (Xl + - X4+) ' (Xl _ - X4")X5 _ X4 + (Xl + - X5 + ) 

(X3~ + 1) X4~ (Xl + - X4 + ) (X4~ + 1) X5~ (Xl + - X5 + ) 



X3" (X4" + 1) (Xl + - X3 + )' § X4" (X5" + 1) (Xl + -X4+)' 

(X5~ + 1) (Xl + + 1) (X2--X4")(X5- + 1)(X2 + -X4~ 

u 9 = —, Z Z^> "10 " 



X5~(Xl + -X5 + ) ' (X4~ + 1) (X2~ -X5") (X2 + -X5 + )' 

(X2~ ~ X5~) (X2 + ~ X5 + ) = (X3~ ~ X5~) (X2 + + 1) (X3 + ~ X5 + ) 

UU (X5- + 1) (X2+ + 1) ' (X2--X5-)(X3 + + 1)(X2+-X5 + )' 

(Xl~ -X5~)(X1 + + 1)X5 + (xi~ + 1) (Xi + + 1) 

u 13 — — f T T\ ) u 14 — Z Z > 

X5 (Xi + - X5 + ) Xl Xi + 

= (Xj~ + 1)X2~ (Xl + - X2 + ) = X2~ (X3~ + 1) (Xl + - X2 + ) 

(Xl"-X2")X2 + ' (X2~ + 1)X3" (Xl + - X3 + )' 

= (X3~ + 1) (X2~ ~ X4~) (X2 + ~ X4 + ) = (X2~ ~ X4~) (X3~ ~ X5~) (X2 + ~ X4 + ) (X3 + ~ X5 + ) 

(X2~-X3")(X4" + 1)(X2 + -X3 + )' 18 (X3" -X4^)(X2" -X5~)(X3 + -X4 + )(X2 + -X5 + )' 

Uig = (X3~ ~ X5~) (X4 + + 1) (X3 + ~ X5 + ) = (Xl~ ~ X5~) (X4 + + 1) X5 + ,g gs 

l!l (X3~ - X4") (X2~ - X5") (X3 + - X4 + ) (X2+ - X5 + ) ' ^ (Xl~ - X4~ ) X4+ (X5 + + 1) ' 

B.2 A T-function parameterization 

Upon expressing the Sij in terms of 20 independent variables q according to the labeling of 
the 20 cross-ratios in the previous subsection, the product Yl^j^^ij) appearing in (5.8) 
becomes 

8 

Y[ r(<%) = r (c 2 - c 3 - c 4 + c 5 + 1) r (c 5 - c 6 - c 7 + c 8 + l) r (c 8 -eg- c w + en + 1) 

i<j 

r (-ci + c 2 + c u - C12 + 1) r (ci 3 ) r (c 6 - c 8 + c 9 + ci 3 ) r (-c 9 - ci 3 - c i4 ) 
r (cu) r (ci + c 9 - en + ci 4 ) r (-ci - c i4 - ci 5 ) r (c 15 ) r (ci - c 2 + c 4 + ci 5 ) 
r (-C4 - ci 5 - c i6 ) r (c i6 ) r (c 4 - c 5 + c 7 + c i6 ) r (-c 7 - c 16 - ci 7 ) r (c i7 ) 
r (c 7 - c 8 + cio + c 17 ) r (-cio - c i7 - cis) r (ci 8 ) r (c w - c u + c i2 + c 18 ) 
r (-C12 - cis - cig) r (cig) r (-c 2 + c 3 + ci 2 + cig) r (-C6 - ci 3 - c 2 o) 

r (-c 3 - cig - c 20 ) r (c 2 o) r (C 3 - C 5 + C 6 + C 20 ) • (B.4) 
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